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Abstract. We investigate the dynamics of a continuous atom laser based on 
the merging of independently formed atomic condensates. In a first attempt to 
understand the dynamics of the system, we consider two independent elongated 
Bose-Einstein condensates which approach each other and focus on intermediate 
inter-trap distances so that a two-mode model is well justified. In the framework 
of a mean-field theory, we discuss the quasi steady-state population of the traps 
as well as the energy distribution of the outcoupled atoms. 
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1. Introduction 

The theory of quantum dissipative systems, and in particular those in which structured 
reservoirs are involved, represents a currently active field of research in a rather broad 
and varied context [I]. The fundamental challenge in problems involving this type 
of reservoir stems from the inapplicability of the Born and Markov approximations, 
normally valid for a reservoir with a smooth density of states, coupled weakly to 
a quantum system with few degrees of freedom. As a result, the so-called pole 
approximation which leads to the elimination of the reservoir degrees of freedom can 
not be adopted. On the other hand, there are no established approaches of general use 
capable of addressing all problems of this type. The situation is further complicated by 
the large variety of and seemingly unrelated physical contexts in which such problems 
may appear. The mathematical structure of the density of states of the reservoir 
and the strength of the coupling to the quantum system are the determining factors 
and source of difficulty. Actually, their combination determines the so-called spectral 
response whose form in each case sets the rules as to which approach may be helpful. 

The theoretical description of what is known as atom laser represents a most 
recent example of this class of problems. In fact, the features that make this problem 
non-Markovian resemble those of an excited atom inside a material with a photonic 
bandgap [2]; even though the underlying physics is rather different. In analogy to 
optical lasers, atom lasers can be obtained by outcoupling atoms from a trapped 
Bose-Einstein condensate (BEC) to free space. The most crucial prerequisite for the 
realization of a continuous atom laser is the pumping mechanism replenishing the 
trapped condensate as atoms are outcoupled from it. A number of sophisticated 
techniques have been proposed to this end [3 , and most of them rely on optical 
pumping between various internal and external atomic states. However, none of these 
techniques has been able to overcome intrinsic losses in the system and achieve laser 
action. Up to date, perhaps the most promising scheme towards the realization of a 
continuous source of condensed atoms was demonstrated by Chikkatur et al. and 
relies on the use of optical tweezers for the transport and the merging of independently 
produced BECs. 

A large number of theoretical models have been used in studies of atom lasers 
which rely mainly on Born-Markov master equations O [6] and Gross-Pitaevskii 
theory [3 [H [9l [TOl [Til HI]- In either case, the models hold only under certain 
operating conditions. For instance, it is well known that for experimentally achievable 
parameters atom lasers may exhibit non-Markovian dynamics |131 114j , which cannot 
be described in the framework of Born-Markov approximations. Besides, inclusion of 
non-Markovian effects in the Gross-Pitaevskii theory is a rather difficult task. Hence, 
investigations of atom-laser dynamics beyond Born and Markov approximations have 
been mainly performed in the framework of a particularly simple model involving a 
single-mode condensate (trap mode) coherently coupled to a continuum of free-space 
modes [El [HI Hg [17l Hg [HI [20l EJ . 

Motivated by the experiments of Chikkatur et al. [4J, in the present work 
we extend these studies to a two-mode scenario. In particular, we consider two 
independent BECs consisting of a large number of bosonic atoms cooled into the 
lowest eigenmode of the corresponding trap. To account for the merging process, the 
two traps are brought together, while atoms are coherently outcoupled from one of 
the BECs only. We focus on an intermediate stage of the merging process where the 
separation of the two traps is so large that, on the one hand a two- mode model can 
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be adopted, while on the other hand a coherent Josephson coupHng is estabhshed 
between the two BECs. Our purpose is to investigate how the presence of the second 
trap mode affects the dynamics of the atom laser and in particular the distribution of 
the outcoupled atoms. 

2. The system 

Our system (see figure [1]) consists of two independently prepared elongated BECs (A 
and B) and let N be the total number of atoms in the system. The two traps are 
initially far apart and each BEC experiences only its local potential V^('B)(r), while 
only the lowest level of each trap (condensate mode) is populated. To allow for the 
merging of the two BECs, the two traps are brought together along one of the tightly 
confining radial directions. Simultaneously, atoms are outcoupled coherently from 
BEC A, by applying external electromagnetic fields. In this section we describe in 
detail the modeling of the system used throughout this work. 

2.1. Double-well potential 

Transport of BECs can be realized using optical tweezers which are produced by 
focused laser beams and offer limited trap volume and depth. Hence, during the 
merging process the two BECs can be brought as close as the traps' beam waist, 
before they start affecting each other. To be consistent with the experimental setup 
for BEC merging [4] as well as related theoretical work [22l [23] , we will assume two 
nearly identical axially symmetric harmonic traps with confining frequencies ujz and 
Trap B is moving towards trap A along the radial direction x (see 
figure [T]) , and the global potential experienced by the trapped atoms can be modeled 
by a time-dependent double-well potential of the form, 

,1 22,1 22 (^\ 

+ -jmujyy + i^muJzZ , (1) 

where m is the atomic mass. According to ([1}, the harmonic potential remains 
unaffected in both y and z directions, while along the merging direction we have a 
double well potential which at any time t, exhibits two minima at x = and x = s{i) 
(see left inset of figure [T]). The distance between the two dips decreases with time and 
at the end of the merging (i.e., a.t t — t^) we have complete overlap. 

Besides the merging time-scale ^m, the details of the motion of trap B are not 
of great importance [H [53]. The crucial point is that the BEC merging must 
be adiabatic so that any kind of excitations in the system are suppressed. To this 
end, first of all the transport of the BECs must take place on a time scale much 
larger than the characteristic time scale of excitations along the merging direction 
i.e., ^ [211[23]. Although this condition can be easily satisfied in a typical 
merging experiment [1] , it does not ensure adiabaticity with respect to the time-scale 
of interatomic interactions [22]. Nevertheless, as long as NUtttm S> 2/iVmodG, where 
Utt is the strength of the interactions and Vmodc is the effective mode volume for each 
trap, it has been shown that only low-lying eigenstates of the Hamiltonian may be 
populated during the merging [22]. As a result, at the end of the process one obtains 
a large single-mode BEC fraction with a unique relative phase. 
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The local potentials centered at x 
equation ([T|) 



and X = s{t), are readily obtained from 



1 



(2) 



In the case of an ideal bosonic gas, the wavefunction of the local ground states |A) and 
|B), corresponi 
like profile i.e. 



|B), corresponding to V^\r) and ^^'"■'(r) respectively, have the well known Gaussian- 



^A(r) = (A|r) = 



1 



^B(r,t) = (B|r) = 



1 
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■ exp 



-s{t)f + U 
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.(3) 



The characteristic harmonic oscillator length Iq is defined as Iq = ^Jh/mLOQ for 
C G {x, y, z}, while l± = Ix = ly As we will see later on, these Gaussian-like profiles 
enable us to obtain analytic expressions for most of the parameters characterizing the 
dynamics of the system. From now on, for the sake of brevity we simply write ipB (r) 
instead of ipB (r, t) . 



2.2. Trapped atoms 

The many-body Hamiltonian describing the dynamics of the trapped atoms is given 

by[nil21[25| 



dY'^l{v) 



-l-V' + vP{v,t) 
2m 



^t(r) + ^ / dr*l(r)*|(r)*t(r)'J't(r), (4) 



where (r) is the annihilation field operator for the trapped atoms with 



#t(r),*l(r') =<5(r-r'), *t(r),^t(r') 



= 0. 



(5) 



The quantity Utt — Anh^att/m measures the strength of the interparticle interaction 
between trapped atoms, while att is the corresponding s-wave scattering length. 

At i = 0, each BEC experiences only its local potential ^('b) as the two traps 
are well separated, and only the lowest level of each trap (condensate mode) is 
populated. We may expand therefore the field operator at i = as ^'t(r, 0) = 
(/?A(r)a(0) -|-i^B(r)&(0), where '/'A(r) and '/'B(r) are the ground-state wavefunctions for 
the traps A and B, respectively [51|211[5S]. The corresponding bosonic annihilation 
operators are denoted by a and h and satisfy the standard commutation relations i.e., 
[di, d]] = (5ij, for di = a and d2 = b. 

To be consistent with the experiment [3] , we will assume that during the merging 
the atomic density profiles follow adiabatically the movement of the traps. The 
condensate wavefunctions start overlapping in space as the two traps approach each 
other and a Josephson-type tunneling is established between the two BECs. If 
the position uncertainty in the ground state of the traps is much smaller than the 
separation of the minima of the global potential V^^*^-'(r, t) i.e., if 

Ix < V2s, (6) 
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the overlap (and thus the Josephson couphng) is small enough so that only the ground 
states of the traps are relevant. In first-order perturbation theory, the corresponding 
local ground-state wavefunctions <^a(b) (r) are orthogonal and describe faithfully BEC 
A and B, at any time < f <C im [23 • Hence, under such conditions we may still 
expand the field operator at times < t ^ tm as [9l [24j [25] 



*t(r,0 =^A(r)a(t) + <^B(r)&(0, (7) 
and the many-body Hamiltonian ^ reduces to the standard two-mode model 

'Ht{t) = hujp^a^a nwBh^h + hJa)h + hJ*P a + hnA_a^a)aa + huBV^y^bb, (8) 
where the coefficients are given by 

u,{t) = i 1 dv^*{v)U^,{v), = ^/ dv\^,{v)\\ Jit) = dr^l{r)L^ipBm) 
with 

Ltit)= -|^V2-fy/«)(r,i). (10) 

In deriving equation ([5]) we have neglected higher-order cross-interaction terms 
involving integrands of the form |iy9i(r)p|((5j(r)p, [ip*{r)ipj{r)]'^, and \(pj{r)\'^ipi{r)ip*{r). 

Condition ^ itself does not justify completely the use of the two-mode model. 
In addition we have to guarantee that the effect of interatomic interactions on the 
ground-state properties of the two wells is small i.e., that hiuixLUyLUzY^^ ^ Nkj [25] . 
For such weakly-interacting bosonic gases, the ground-state wavefunctions are well 
approximated by equations ([3]) . The Gaussian profile of the wavefunctions enables us 
to evaluate analytically all of the integrals ^ obtaining for the coefficients entering 
the Hamiltonian Tit 

UJj = LUo, Kj = K, J = J*, (11) 



where 



UJo{t) = UJz 



' ' ' ■ ^Erfc(^) " 



2 A A Ay? 
11 r; ' 
2 ^ A " AV^, 



(12) 



.m^u;Al + \--^]e-^\ (13) 



hat, 



ATOV27rZ? 



(14) 



while Erfc(ry) is the complementary error function and we have introduced the 
dimensionless quantities r]{t) = s{t)/{2lx) and A — ujzjuJx- Condition h{u:x^yUJzYl'^ ^ 
iV^j thus yields the following upper bound on the total number of atoms we may 
consider 

7V« Ai/3\/2^-^. (15) 
att 

For the reasons we discussed earlier in this section, throughout our simulations 
we focus on inter-trap distances > 1.5. According to (fT3|) . J{t) is practically zero for 
large 77 and increases in absolute value, as we bring the traps closer i.e., for decreasing 
r\ (see right inset of figure [T|). On the contrary, the ground-state frequency tJo(i) does 
not vary appreciably in the same regime of inter-trap distances and thus throughout 
our simulations we may safely assume that lOq remains practically constant i.e., 

^o{t) «^ + ^) • (16) 
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Finally, before we proceed further, it is worth recalling here that the two 
lowest eigenstates of the global double-well potential can be well approximated as 
the symmetric and antisymmetric combinations of the local eigenstates i.e., |±) = 
(|A) ± |B))/\/2, with eigenfrequencies uj± = Wo ± J- Hence, Hamiltonian ([5]) can be 
also expressed in terms of the corresponding symmetric and antisymmetric bosonic 
operators d± = (a ± 6)/V2. (e.g., see [25]). Actually, in some cases one might get 
further detailed insights into the dynamics of the system if this is viewed in the basis 
of the global states. 

2.3. Outcoupled atoms 

Let us assume that atoms are coherently coupled out of the BEC A only. Neglecting 
collisions between trapped and free atoms, the many-body Hamiltonian for the free 
atoms is simply of the form [9] 

Ht^ J dr¥,{r)Lti>i{r), (17) 

where 

Lf^~^V + Vfir). (18) 
Zm 

In general, the potential Vf(r) experienced by the free atoms depends on the 
particular setup under consideration. Throughout this work we consider an atomic 
waveguide for the outcoupled atoms [241 126] , resulting in an effective one-dimensional 
atom laser propagating along the weak confining axis of the waveguide (see figure [1]). 
For instance, such a guided atom laser has been demonstrated recently by Guerin et al. 
|27| and offers many advantages over the conventional outcoupling schemes. Formally 
speaking, the strong transverse confinement allows us to assume that the transverse 
dynamics of the free atoms adiabatically follow the slowly varying transverse potential 
of the optical guide Vf (rj^) [27]. For the sake of simplicity, throughout this work we 
assume that the transverse guide potential is nearly the same with the transverse 
potential of trap A i.e., V{{r±) ~ V^'^-'(rj^). In the absence of gravitational or other 
forces (as in the experimental setup [27 ) , the longitudinal component of the potential 
is V[{z) = 0. Thus, the field operator for the free atoms can be expanded as 

*f (r, t) = (^A(r±) J2 Xk{z)ck{t), (19) 

k 

where Ck is the annihilation operator of free atoms with momentum hk and obeys 
the usual bosonic commutation relations [cfc,c^,] = Sk,k'- The wavefunction (/?A(r±) is 
the ground-state wavefunction of the local transverse potential V^'^''(rj_), with the 
normalization J dr±\(pA{r±)\'^ = 1, so that the linear atomic density piuizjt) = 
J dr±\^i{r,t)\'^ — J2k qXk^q'^k'iq- The longitudinal wavefunction Xk{z) is readily 
obtained as a solution of the time-independent Schrodinger equation for a free 
atom (i.e., for V{{z) = 0). Thus for a free atom with momentum hk we have 
Xkiz) = e^'^^/V^TT, and frequency 
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As we will see later on, this quadratic dependence of LUk on k is responsible for a 
number of mathematical difhculties arising in the context of atom lasers [13]. Using 
expansion (|19p and the orthonormality condition for Xk{z), Ti-i reads 

nf^hJ2{^k + '^)cick. (21) 



2.4- Output coupling 

We consider an output coupling by application of external electromagnetic fields which 
induce an atomic transition from the internal state (\t)) of the trapped atoms to an 
untrapped state |/). In the rotating- wave approximation, the many-body interaction 
Hamiltonian is of the form 9J 

Vit) = hj dr#t(r) ^A(r,t) ^-t^ + H.c (22) 

where A(r, t) is the coupling between trapped and untrapped atomic states. Using the 
expansions ([7]) and (fT9)) . we obtain 

V{t) = -^J2ci f drv/AM)^l(r^)e-''=^ U(r)a(t) + ^B(r)6(t)l + H.c (23) 

In general the form of A(r,f) depends on the particular outcoupling mechanism 
under consideration. Typical mechanisms may involve one-photon radio-frequency 
transition or indirect two-photon stimulated Raman transition [II [3 [9l [Til [El [26] . 
Note that in the following we neglect the momentum kick experienced by the atoms 
as well as the spatial dependence of A(r,t), obtaining 

V(t) = hY^g{k,t) (ac^ a^ck) + fie""' ^5(fe,t) (bcl + b^Ck) ,(24) 

k k 

where 

gik,t) ^ ^^^)e-^'^V\ (25) 

According to equation ([M]) , the interaction consists of two terms despite the fact 
that the outcoupling mechanism is applied on BEC A only. More precisely, the first 
term of V{t) refers to BEC A and is similar to the expression used by many authors in 
the context of the standard single-mode model for the atom laser [T31 [HI [HI [T71 [TSl 
[Tni[2ni[ll]- However, due to the presence of the second BEC, in our model we have 
obtained one more term which is proportional to the overlap of the two ground-state 
wavefunctions, that is e~'' . In that respect, equation ([24)) is a generalization of the 
single-mode outcoupling Hamiltonian [ni[TSl[TCl[T71[TSl[Tni[2ni[2T], toa two- mode 
scenario. 

Let us now estimate the spectral response of the atomic continuum for the 
particular outcoupling mechanism under consideration. The density of states which 
are available to a free atom can be determined by the dispersion relation ([20)) as follows 

dk 



p{uj) = 



dio 



-,/^e(.), (26) 



where 6(u;) is the usual step function. Note the divergence of the atomic density of 
states at the edge frequency Wo = 0, which is a characteristic property of the one- 
dimensional model under consideration. Taking advantage of the symmetrical shape 
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of the coupling and the even parity of uj^, we may reduce the fc-space only to the A; > 
sub-space. The spectral response of the continuum is then of the form 



where g{ijJ^ t) is readily obtained from g(k, t) using the atomic dispersion relation (j20|) . 

At this point we have completed the presentation of our model and the underlying 
approximations. In closing, let us summarize the main results by rewriting the 
complete form of the Hamiltonian under consideration in a frame rotating at lo±_ 

n = —{a'' a + b^'b) + w^c^Cfc + hK{a)a^aa + S^S^SS) + hj{a)b + V^a) 

k 

+hJ29ik, t){acl + a^Ck) + ^ g{k, t){bcl + b'^cu). (28) 



3. Heisenberg equations of motion 

Given the total Hamiltonian (j28p one may proceed to derive Heisenberg equations of 
motion for the operators of interest. In the Heisenberg picture, the evolution of the 
expectation value of an arbitrary operator A is governed by 

Thus, for the operators pertaining to the two traps and the continuum, we obtain 



d(a) 
dt 



= ~i—{a) -2iK{a)aa) -iJ{b) -2i I dkg{k,t){ck) , (29) 







dt 

d{ck 



- i'^ib) - 2m{Pbb) - U{a) - 2ie-''" dkgik,t){ck) ,m 

- itokick) - igik,t){a) ~ igik,t)e-'^\b). (31) 



dt 

We may now distinguish between two cases. 

In the absence of interatomic interactions (i.e., for k — 0) the Hamiltonian 
(|^ becomes bilinear. As a result, the above set of equations is closed and all of 
the initial statistical properties of the system are preserved in time. For instance, 
if the BECs are initially in coherent states, we have (a^(t)a(t)) — (a^(t))(a(t)), 
(St(t)6(t)) = {b^{t))(b{t)) and {cl{t)cjit)) = {c\{t)){cj{t)), for all t > (see also 
[T51 [TBI [T71 [5T]). In other words, the bilinear form of the Hamiltonian preserves the 
initial coherence in time, so that at any instant t we can decorrelate exactly any 
higher-order correlation function in terms of (a), (6), and {ckj- 

In the presence of interatomic interactions (i.e., for k 7^ 0) the Hamiltonian (|28p 
involves fourth-order terms, and thus we have the appearance of third-order correlation 
functions in the right-hand side of equations (P^ - ([?T|) . This set of equations is 
no longer closed, while consideration of differential equations for the third-order 
correlation functions leads to the appearance of terms of even higher order and so 
on. In general, there are no exact remedies for such mathematical problems, but an 
approximate solution can be always obtained by decorrelating higher-order correlation 
functions into products of lower ones. 
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In the present work we decorrelate the third-order correlation functions appearing 
on the right-hand side of equations (|29p - (PT|) as follows: {a^aa) w {a^){a){a) and 
{b''bb) « {b^){b){b). Hence, equations jlUl-lISl]) read 
d(a) .uj^ 



— 1- 



-(a) -2iK|(a)na) -iJ(6) -2i / dkg{k,t){ck) , (32) 







dt 2 

^ = -i^{b) -2iK\(b)\^(b) -U {a) -2ie-^' dkg{kMcim 

^ = -\0Jk{ck)-^9{kMa)-ig{Kt)e-^\b). (34) 

One way to solve such a set of coupled differential equations is by means of the 
Laplace transform method. To this end, however, one has to be able to perform all 
of the integrations over the continuum as well as the inverse Laplace transforms at 
the end. Both of these tasks are more or less straightforward in the case of smooth 
continua for which the Born and Markov approximations are applicable. In the present 
context, however, none of the aforementioned approximations is valid. Indeed, as we 
discussed earlier, the quadratic atomic dispersion relation is associated with a density 
of atomic states which diverges for small frequencies (see equation [26|) . This behavior 
is also reflected in the spectral response ((27|) and implies that the continuum under 
consideration does not vary slowly for all frequencies. 

Structured continua which invalidate both Born and Markov approximations 
emerge in different areas of physics and have attracted considerable interest over 
the last few years [1]. To address fundamental mathematical difficulties associated 
with these continua a number of new theoretical techniques have been developed 
[U [ini [28l [29l [30] . Here, to deal with the structured continuum at hand, we follow a 
discretization approach developed in the context of photonic band-gap continua |31j . 
Briefly, we substitute the continuum for frequencies within a range around lOq (i-e., 
for < a; < Wup), by a number (say M) of discrete modes, while the rest of the 
atom-mode density is treated perturbatively since it is far from resonance. Discussion 
on the choice of t^up and the number of discrete modes can be found in Refs. |31j . 
This approach has been also applied in the context of atom lasers |5T] , and is capable 
of providing not only the evolution of the number of atoms in the condensates, but 
also the distribution of the outcoupled atoms in frequency domain, irrespective of the 
strength of the outcoupling and the form of the spectral response. 

In general, a continuum can be discretized in many different ways (see for instance 
[21j ) and in the present work we have chosen a uniform discretization scheme. In 
particular, we choose the frequencies of the discrete modes to be Ljj = je, where 
the mode spacing e is determined by the upper-limit condition of the discretization, 
namely Wup — Me. The corresponding coupling for the j mode, is determined by the 
spectral response (|27l) as follows 

g] = D{cj,)e. (35) 

Hence, working similarly to [H], equations ) read after the discretization 



M 



dt 



^^^^ - ■■ - S) {a) --2i^\{a)\^a) -i{J - Se-^->'){b) -iY.~g,{c,), (36) 

i=i 

M 

_ Se-'^') {b) - 2i^|(&)p(&) - i(J - Se-^>'){a) - ie-'^ ^.g,(c,)(37) 



dt 
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= - i^^j(cj) - igj(a) - iffj-e (b), (38) 



OO 



/ 2^du. (39) 



where 



4. Simulations 

Throughout our simulations we have considered ^'^Na BECs with m = 3.818 x 10~^^ 
Kgr and att = 2.75 x 10~^ m, which are formed independently in identical harmonic 
traps with longitudinal oscillation frequency cj^ — 200 sec~^ and ratio A = 0.4. We 
assume that the BECs A and B are initially prepared in coherent states \a) and |/?), 
respectively. Equations ([36 |) -([38 |) are thus solved with initial conditions 

{a{0))^a^^/Nm, (6(0))-/3- V^iV/3(0)e'*W, (c,(0))=0, (40) 

where 0(0) is the initial relative phase between the two BECs. Accordingly, the 
initial number of condensed atoms in the traps A and B are given by iVA(O) = 
|(a(0))|2 = 7V5(0) and iVB(O) = |(6(0))P = iV/3(0). At any time t > we have 
a{t) + /3(i) + 7(t) = 1 so that A^a(0 + ^b(0 + Nc{t) = N , where Nc{t) and 7(t), 
respectively are the population and amplitude of the continuum. Finally, for the sake 
of simplicity and without introducing any significant errors, the applied outcoupling 
pulse A(i) is modeled as rectangular lasting from t = to t = r. 

Most of the work on the non-Markovian aspects of atom-laser outcoupling has 
been performed in the framework of an ideal gas (i.e., for hi — Q) [T31 [121 [HI [IZl [ISl 
[T9l [20l [2T] . Hence, for the sake of comparison, in this section we focus mainly on 
the analysis of results obtained by simulations in non-interacting systems. The case 
of weakly-interacting gases (k ^ 0) deserves a thorough investigation and as such 
will be discussed in detail elsewhere. At the end of this section, however, we briefly 
summarize some of the main features of the weakly-interacting systems we have found 
in our simulations. 



^■l. Ideal bosonic gas 

As we mentioned earlier, in the absence of interatomic interactions no decorrelation 
approximations are necessary for the derivation of a closed set of equations for the 
expectation values of the operators {a{t)), (6(t)) and {cj{t)). The evolution of the 
system is obtained by propagating equations p6)) - (|38)) with k = 0. As both BECs are 
assumed initially prepared in coherent states, all of the statistical properties of the 
system at any time t can be expressed exactly in terms of (a(t)), {b{t)) and {cj{t)). 



4.1.1. Weak outcoupling — Markovian dynamics. As depicted in figure [21 for weak 
outcoupling strengths (i.e., for A < 5 x 10^ sec^^), the dynamics of the system are 
mainly Markovian. More precisely, we have population exchange between the two 
traps, but the oscillations are exponentially damped as atoms are irreversibly coupled 
out of the traps. As we reduce the distance rj between the traps, the Josephson 
coupling J increases and the oscillations become faster (e.g., compare figures [2][b) and 
[21(c)). On the contrary, for constant intertrap distance, the oscillations decay faster as 
we increase the outcoupling rate A (compare figures [2Kb) and[2Kd)). 
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In general, the distribution of the outcoupled atoms (see figure [3]) exhibits the 
characteristic doublet which, however, is expected to be asymmetric mainly due to the 
unconventional density of atomic states (pS)) . In particular, the origin of the doublet 
is well-described in terms of the global states |±) with eigenfrequencies 0)2/2 ± J. The 
intertrap coupling splits the previously degenerate local states |A(B)) into a doublet of 
global states and thus the outcoupled atoms emerge as distinct peaks separated by 2J . 
The outcoupling rate is larger for the global state with frequency closer to the edge 
Uc = 0, as the density of availlable atomic states, and thus the spectral response, scales 
as l/V^- Decreasing the distance between the two traps, the Josephson coupling 
increases monotonically for 77 < 2.0, and thus the distance between the peaks also 
increases (compare figures [S^b) and [Sl^c)). On the other hand, as depicted in the 
inset of figure[Tl for 77 > 2.0 the Josephson coupling is not strong enough to produce a 
noticeable splitting and to give rise to a clear doublet in the distribution of outcoupled 
atoms (see figure Elja)). Finally, as we increase the outcoupling rate A for constant 77, 
the peaks become broader while their position remains practically unchanged (compare 
figures Hb) andUd)) . 

4-1-2. Strong outcoupling — Non-markovian dynamics. In the strong-outcoupling 
regime (i.e., for A > 5 x 10^ sec~^) the trapped populations begin exhibiting non- 
Mar kovian dynamics. In general, the evolution of the system is governed by two 
different processes namely, the exchange of population between the BECs and the 
exchange of population between the BECs and the continuum. Note that exchange 
of population between any discrete feature and a continuum is a signature of the 
non-Markovian nature of the problem under consideration. 

In figure |4] we present the evolution of the trapped populations as functions of 
time, for a given intertrap distance and increasing outcoupling rate. Clearly, as far as 
trap B is concerned, we can identify an initial transient regime where the main part 
of the population is lost. After this initial stage, dissipation is temporarily turned 
off and trap B gets atoms from trap A (slight oscillations). This weak oscillatory 
population exchange between the two traps persists even for larger times, but the 
population of trap B is gradually transferred into the continuum in an irreversible, 
almost exponential, way. Thus, irrespective of the strength of the outcoupling rate, 
trap B is always empty in the long-time limit (e.g., see figures |4Ka,b)). 

Moreover, in figure [4] we see that, besides the weak oscillations, the population 
of trap A exhibits fast oscillations which become more pronounced and faster with 
increasing outcoupling strength. These oscillations are reflected only in the population 
of the continuum which is not shown here. For relatively weak outcoupling rates not 
only trap B, but also trap A is empty in the long-time limit (e.g., see figures |4Ka,b)). 
On the contrary, for stronger outcoupling rates the system reaches a steady state 
pertaining to a practically empty trap B and a partly depleted BEC A (e.g., see figures 
|4l[d)). The formation of such a bound state has been demonstrated experimentally 
[14j and involves atoms in a superposition of two states namely, the trapped and 
the untrapped state. The Born and Markov approximations are valid only if such 
superpositions decay on a time scale much shorter than the time scale of interest in 
this work. It is also worth noting that, according to our simulations (e.g., see figureH]), 
the bound state involves only trap A and not trap B. This is perhaps due to the weak 
outcoupling rate for BEC B which is proportional to the overlap e"** between the 
two BEC wavefunctions. Hence, BEC B is only weakly coupled to the continuum and 
any non-Markovian effects, such as the formation of a bound state, are suppressed. 



Non-markovian dynamics in atom-laser outcoupling from a double-well BEC 12 

Finally, although limitations in the validity of the two-mode model does not allow 
us to consider larger values of A, from figure Elja) it is obvious that the population 
trapping increases as we increase the outcoupling rate. 

Let us discuss now the effect of the intertrap distance rj on the evolution of the 
trapped populations. As depicted in figure[6](a), for large intertrap distances the weak 
oscillations are absent and the population of trap B remains practically constant as 
the overlap between the BEC wavefunctions is negligible. As we bring the traps closer 
(see figures [6] (b-d)), the population of trap B starts evolving in time with an initial 
transient regime followed by an irreversible decay. On the other hand, for the typical 
values of 77 allowed by our two-mode approximation, there seems to be no significant 
effect on the evolution of the population of trap A. In particular, we have a continuous 
exchange of population between trap A and the continuum until the formation of a 
bound state in the long-time limit. Although the steady-state population of trap A 
does not change considerably over the regime of intertrap distances rj we can cover, a 
slight reduction is noticeable in figure[5Jb) as we reduce rj. Moreover, there is a regime 
of distances around rjc ~ 2.14, where the system seems to have no steady state, in the 
conventional sense of the word, as it is beating between the two localized condensate 
modes even for the long time scale used in figure [3 That is why we do not give any 
values for {a) {T)a{T)) in the neighborhood of rjc. All of this behavior can be easily 
understood if the system is viewed in the basis of the symmetric and antisymmetric 
global states |±). 

First of all let us briefly discuss the behavior of the global states |±) with 
eigenenergies lo± — ujz/'i ± J, as we approach the two traps. As depicted in the 
inset of figure [1] the Josephson coupling and thus the separation of the global states, 
does not vary monotonically with rj. More precisely, as we reduce rj the splitting of 
the states |±) increases for 77 > ?7inin, where Tymin ~ 2.3. In this regime of intertrap 
distances, J is negative and thus lj- > lu^ i.e., the symmetric state moves to the 
left and the antisymmetric state to the right of the characteristic frequency uJz/2. 
This relative movement is inverted for rjc < rj < rj^nin and the two global states 
begin approaching each other. They become basically resonant ai rj — rjc where J 
changes sign. Subsequently, as we further reduce rj, the splitting of the states |±) 
increases again, but this time the global states move in opposite directions as J > 
and w+ > LO-. In particular, the symmetric state moves away from the edge and 
the population trapping associated with it decreases, whereas the antisymmetric state 
moves towards the edge where it is more protected against dissipation. In any case, it 
is worth noting that the Josephson coupling attains significant values only for 77 < 2.0. 
Thus, for 2.0 < rj < 3.0 the two condensate modes are practically resonant, as the 
splitting is negligible. Hence, the system is beating between the two condensate modes 
and no steady-state is found on the time scale of figured 

In view of the above discussion, it is also easy to understand the reduction of the 
steady-state population of trap A as we bring the traps closer. For the particular 
initial conditions under consideration, we obtain (d!^(0)d+(0)) = \a + f3\^/2 and 
(dL(0)d_(0)) = |a — /?P/2. Throughout this work we have focused on in-phase 
BECs only i.e., = 0. As a result, the main part of the population initially 
occupies the symmetric global state and thus the behavior of the corresponding steady- 
state population with decreasing rj also determines the behavior of the steady-state 
population in the local state |A). This may not be the case if we choose different 
initial conditions, but the effect of the phase difference on the system's dynamics 
will be investigated in detail elsewhere. 
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We turn now to the discussion of the distribution of the outcoupled atoms in the 
non-Markovian regime. Unfortunately, for the spectral response (P7)l the derivation 
of analytic expressions for the atomic distribution is a rather difficult task [Tn|. In 
fact, analytic results can be obtained only in some special cases e.g., in the limit of 
broad-band output coupling, by means of the Laplace transform method [16]. The 
discretization approach, however, is capable of providing us with the distribution of the 
outcoupled atoms at any time. In figures [7| and [8] we present such typical distributions 
at the end of the outcouphng pulse, i.e. at times t — t, with the pulse duration chosen 
sufficiently large to ensure that the distributions do not vary significantly with time. 

For a better interpretation of these results, it is worth keeping in mind that the 
system under consideration involves two condensate modes which decay into the same 
atomic continuum. Moreover, there are two different outcoupling channels for each 
mode. More precisely, atoms can be coupled out of BEC A either directly or via BEC 
B and vice- versa; albeit at different rates. As a result we expect quantum interference 
effects which in addition to the non-Markovian nature of the dynamics may give rise 
to unconventional distributions of the outcoupled atoms. 

When the traps are far apart (see figure [TlJ a)), we are essentially in the single- 
trap case where atoms are coherently outcoupled from BEC A only. In this case any 
quantum interference effects are absent and the distribution of the outcoupled atoms 
exhibits a well-known profile previously discussed in Ref. [H]. More precisely, we 
have a peak around the condensate-mode frequency uJz/2, and a peak around the 
edge frequency We = where the density of atomic states diverges (see equation [26]) . 
The former peak is shifted towards higher frequencies due to the coupling between the 
condensate mode and the continuum. This effect has also been noted by other authors 
(e.g., see [20]) and becomes more pronounced as we increase the outcouphng rate. 
Moreover, the time-dependence of the presented distributions gives rise to oscillations 
which become faster as we increase the pulse duration. In the Markovian regime, 
the atomic distributions do not exhibit any such oscillations (see figure |3]) as all of 
the time-dependent terms entering the distributions become practically negligible for 
times At ^ 1. In the non-Markovian regime, however, we have a continuous exchange 
of atoms between the BECs and the continuum even for larger times as well as the 
formation of a non-decaying bound mode in the long-time limit. Hence, there exist 
time-dependent terms associated with the formed bound mode which persist even in 
the limit t —^ oo and their effect on the atomic distribution is evident, unless one 
performs a time average over the period of the oscillations (see also discussion in 

[S1IS2). 

As the traps approach each other (see figure [Tjjb)), the overlap between the two 
BECs increases and a narrow peak appears next to the main peak. This narrow peak 
is clearly associated with a rather weak outcoupling and thus can be attributed to 
atoms originated from trap B. It could be said therefore that for intermediate intertrap 
distances the distribution of the outcoupled atoms is basically a superposition of the 
distributions for each individual trap. As we reduce further the intertrap distance rj, 
quantum interference effects start becoming significant, and a clear dip appears next 
to the narrow peak (see figures [7l^c,d)). 

As depicted in figure [DJa) , the dip is also present for relatively weak outcoupling 
strengths provided the two traps are close enough. Actually, it was also present in 
the Markovian regime (see figure [3]), but it was not so clear as the two peaks were far 
apart. It is obvious therefore that the dip is a clear evidence of destructive interference 
between the various outcoupling channels of the system. Although the outcoupling 
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itself is not sufficient to cause the dip, in figure [H we see that the dip becomes broader 
with increasing outcouphng rates, while one may also notice the emergence of the 
narrow peak at frequencies around uJz/'i in figures [5]Jb)-(d). 

Due to the lack of analytic expressions, it is not clear whether the observed dip 
at LOj = is a perfect dark line i.e, whether {<^j{T)cj{T)) = at uj^. To resolve 
this issue we have obtained analytic expressions for the distributions in the broad- 
band limit of the output coupling where -D(w) ~ h.{t)/\/io [HI [16]. In this limit, the 
observed profiles are similar to the ones presented here and we can verify that the 
atomic distribution at ujj — lo^ vanishes only in the limit t ^ oo. Otherwise, for 
M:^ I the dip corresponds to a very low probability (~ 10^^) for outcoupled atoms 
with frequency cJd- 

4-2. Weakly-interacting bosonic gas 

Considering a harmonic trap with longitudinal oscillation frequency oj^ — 200sec~^ 
and ratio A = 0.4, condition yields N <C 2.5 x lO'^. In other words, our model is 
valid for small BECs consisting of a few hundred of atoms. As we discussed in scction[3l 
in the presence of interactions we can obtain a closed set of equations of motion for the 
operators {a{t)), {b{t)) and (cj(i)), only by applying a decorrelation approximation. In 
view of this decorrelation, interatomic interactions enter the equations of motion for 
{a(t)) and {b{t)) as time-dependent shifts proportional to the corresponding trapped 
populations TVaIO = l(a(0>P and Nb = (see equations [36l and [37]) . As a 

result, the frequencies of the condensate modes A and B fluctuate in time and become 
off-resonant. 

The evolution of the weakly interacting system is governed by three distinct 
physical processes. More precisely, apart from the Josephson and the output coupling 
which were also present in the interaction-free model, we also have the repulsive 
coUisional interactions. It is reasonable therefore to define the ratios Nc^^ — J/k and 
ivi^^ = VH/k which quantify the effect of interatomic interactions relative to tunneling 
and outcoupling effects, respectively. How strongly the inclusion of interactions affects 
the results obtained in the framework of the interaction-free model depends on the 
ratios N/Nc^'' and N/Nc\ where N is the total number of atoms in the system. 

The case of weakly-interacting gases is of particular interest but it cannot be 
covered in the present work as there are many aspects which need to be thoroughly 
investigated. For instance, it is already known that for an isolated double-well BEC 
(i.e., in the absence of losses and outcoupling) one may define two extreme regimes 
of dynamics [34 namely the Josephson regime [35| and the self-trapping regime |34| 
. The detailed analysis of the previous section does not involve self-trapping at all as 
this phenomenon occurs only for Hamiltonians involving interactions. Here, we would 
like to briefly highlight only some of the features of the weakly-interacting model we 
have found in our simulations. The detailed presentation and discussion of the results 
will be the subject of a forthcoming work. 

In the strong outcoupling regime (i.e., for A^c > A^), inclusion of interactions 
affects the evolution of the trapped populations only quantitatively. The most 
important features of the weakly-interacting system in this regime seem to be the 
destruction of the bound mode and the disappearance of the dark spectral line 
discussed in the context of the interaction-free model. On the other hand, for Nc < N 
interatomic interactions dominate over the output coupling and thus we have strong 
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modifications in the evolution of the populations as well as the distribution of the 
outcoupled atoms. 

5. Summary and outlook 

We have investigated the non-Markovian aspects of atom-laser outcoupling from a 
double- well BEC. Our two- mode trapped condensate model has been motivated by 
recent experiments on the merging of independently formed BECs and the first 
realization of a guided quasicontinuous atom laser [27] . For the sake of comparison 
with earlier work, relying on a single-mode trapped condensate, we have focused on 
an interaction-free model. In particular, we have studied how the presence of the 
second BEC (BEC B) affects the evolution of the trapped populations as well as 
the distribution of the outcoupled atoms. Although the outcoupling mechanism is 
basically applied only to BEC A, atoms are also weakly outcoupled from BEC B due 
to the overlap between the BEC wavefunctions. 

In the case of weak outcoupling rates the dynamics of the system are purely 
Markovian and thus particularly simple. More precisely, the system oscillates between 
the two condensate modes while atoms are coherently outcoupled into the continuum. 
The distribution of the outcoupled atoms exhibits the characteristic asymmetric 
doublet due to the unconventinal density of atomic states and the different outcoupling 
rates experienced by the BEC atoms. 

The situation is substantially different in the strong-outcoupling regime, where the 
system exhibits non-Markovian dynamics. In particular, we have population exchange 
between the two BECs as well as between the BECs and the continuum. In the 
latter process, BEC B seems to participate passively as its population is gradually 
transferred to the continuum in an irreversible way. On the contrary, BEC A keeps 
exchanging population with the continuum even for larger times while for sufficiently 
strong outcoupling rates it is found only partially depleted in the long-time limit. The 
formation of such a bound state has also been predicted in the context of single-mode 
models. However, in the two-mode model under consideration the long-time-limit 
behavior of the system seems to depend on several parameters such as the outcoupling 
rate, the intertrap distance and the phase difference between the two BECs. 

The non-Markovian nature of the dynamics and the presence of the second BEC 
are mostly apparent in the distribution of the outcoupled atoms which exhibits two 
peaks. On the one hand, there is a broad peak stemming from atoms directly 
outcoupled from BEC A while on the other hand, atoms outcoupled from BEC 
B give rise to a rather narrow peak. Most importantly, as a result of destructive 
quantum interference between various outcoupling channels in the system, the atomic 
distribution may also exhibit a dark spectral line. It is worth noting that similar 
quantum interference phenomena have been discussed in the context of optical 
systems. Nevertheless, to the best of our knowledge the appearance of dark lines 
in the spectrum of atom lasers has not been discussed in the literature so far. 

In general, the two-mode model considered here allows for the direct inclusion of 
collisional interactions between trapped atoms. Nevertheless, due to space limitations, 
throughout this work we have focused on an interaction-free model only, and have 
assumed that the two BECs are initially in phase. Effects of interactions and the role 
of the phase difference will be presented in detail elsewhere. 
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Figure captions 




Figure 1. Schematic representation of the system under consideration. Left 
inset: The double- well potential experienced by the trapped atoms along the x- 
direction. Right inset: The Josephson coupling as a function of the dimensionless 
intertrap distance rj. 
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Figure 2. Markovian regime. Evolution of the normalized trap populations as 
a function of time for various inter-trap distances and outcoupling strengths: (a) 
A = 10^ sec-2, ri = 2.0; (b) A = 10^ sec'^, r) = 1.7; (c) A = 10^ sec'^, = 1.5; 
(d) A = 2 X 10^ sec~2, rj = 1.7. System parameters: u)^ = 200 sec~^, A = 0.4. 
Initial conditions: a(0) = 0.7, /3(0) = 0.3. Discretization parameters: M = 1500, 
<.^up = 300 scc^-"^. 
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Figure 3. Markovian regime. Distribution of the outcoupled atoms at r = 10 sec 
for various inter-trap distances and outcoupling strengths. The distributions (a)- 
(d) are in one-to-one correspondence with figures[2la)-(d), respectively. Note that 
for the sake of illustration the regime on the left-hand side of the vertical dashed 
line has been magnified ten times. 
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Figure 4. Non-markovian regime. Evolution of the normalized trap populations 
as a function of time for various outcoupling strengths: (a) A = 5 X 10^ sec~^; (b) 
A = 10^ sec~^; (c) A = 2 X 10^ sec~-^; (d) A = 4 X lO"* sec~^. System parameters: 
LJz = 200 sec-l, A = 0.4, rj = 1.7. Initial conditions: a(0) = 0.7, /3(0) = 0.3. 
Discretization parameters: M = 1500, uJup = 300 sec~^. 
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Figure 5. Non-markovian regime. Typical behavior of the steady-state 
population of trap A for varying outcoupling rate (a) and dimensionless intertrap 
distance (b). The depicted values are estimated numerically for a sufficiently dense 
discretization and pulse duration r = 40 sec. In the neighborhood of »? f» 2.14, we 
have not found a steady-state for trap A, although we have let the system evolve 
for times up to r = 50 sec. System parameters: uiz = 200 sec~^, A = 0.4, t] = 1.7 
for plot (a) and A = 4 x 10^ sec~^ for plot (b). Initial conditions: 5(0) = 0.7, 
/3(0) = 0.3. Discretization parameters: M = 3000, w^p = 300 sec~^. 




Figure 6. Non-markovian regime. Evolution of the normalized trap populations 
as a function of time for various intcr-trap distances: (a) rj = 4.0; (b) r; = 2.0; 
(c) r; = 1.8; (d) rj = 1.6. System parameters: lUz = 200 scc^-', A = 0.4, 
A = 4 X 10^ sec~^. Initial conditions: a{0) = 0.7, /3(0) = 0.3. Discretization 
parameters: M = 1500, uiup = 300 sec~^. 
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Figure 7. Non-markovian regime. Distribution of the outcoupled atoms at 
T = 10 sec for A = 2 X 10^ sec"'^ and various inter-trap distances: (a) r} = 4.0; (b) 
rj = 2.0; (c) r] = 1.8; (d) r] = 1.6. System parameters: u)z = 200 sec~^, A = 0.4, 
A = 2 X 10^ sec-2. Initial conditions: «(0) = 0.7, fifSi) = 0.3. Discretization 
parameters: M = 1500, ojup = 300 sec~^. 




Figure 8. Non-markovian regime. As in figure [7] for r) = 1.7 and various 
outcoupling strengths: (a) A = 5 X 10^ sec~^; (b) A = 10"^ sec~^; (c) 
A = 2 X 10^ sec-2; (d) A = 4 X 10^ sec-^. 




Figure 9. Non-markovian regime. A closeup of the atomic distribution around 
the dip for the parameters of figure \E\ and various outcoupUng strengths: A = 
5 X 10^ sec~2 (solid Hne), A = 10^ sec"^ (dotted line), A = 2 X 10^ sec"^ (dashed 
Une), and A = 4 x 10^ sec~^ (dot-dashed hne). 



